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Abstract 

cn 
o 

Lines in the energy spectrum of gamma rays are a fascinating experimental signal, 
which are often considered "smoking gun" evidence of dark matter annihilation. The 
current generation of gamma ray observatories are currently closing in on parameter 
^ space of great interest in the context of dark matter which is a thermal relic. We 

consider theories in which the dark matter's primary connection to the Standard Model 
is via the top quark, realizing strong gamma ray lines consistent with a thermal relic 
through the forbidden channel mechanism proposed in the Higgs in Space Model. We 
consider realistic UV-completions of the Higgs in Space and related theories, and show 
that a rich structure of observable gamma ray lines is consistent with a thermal relic 
as well as constraints from dark matter searches and the LHC. Particular attention 
is paid to the one loop contributions to the continuum gamma rays, which can easily 
swamp the line signals in some cases, and have been largely overlooked in previous 
literature. 



1 Introduction 



The existence of nonbaryonic dark matter is well established, but there is still no clue as 
to its particle identity, and no indication as to what interactions beyond gravitational it 
experiences. In order to place dark matter in context within the Standard Model or its 
extensions, observation of dark matter interacting through something other than gravity 
will be required. There is currently a rich and diverse experimental program searching for 
dark matter through its direct scattering with nuclei, produced at colliders, or indirectly 
through its annihilation products. 

This last pillar of the searches for dark matter is particularly interesting. Observation of 
dark matter annihilation would be a clear indication that dark matter has non-gravitational 
interactions, and (depending on the observations) would help to establish it as a thermal 
relic. In fact, observations of dwarf spheroidals by the Fermi LAT currently exclude s-wave 
annihilating thermal relic dark matter with mass below about 50 GeV [lj[2]- 

An important element in the indirect search for dark matter is the search for lines in 
the energy spectrum of gamma rays. Such features occur when dark matter annihilates 
into a two-particle final state, with one of the particles a photon. Since dark matter must 
be electrically neutral (to good approximation [3]), such a signal is typically expected to 
be rather smaller than the expected annihilation cross section of a thermal relic, (crv) ~ 
3 x 1CT 26 cm 2 /s (with mild dependence on the dark matter mass jij). Nonetheless, such a 
gamma ray line is a very distinctive signal, unlikely to be faked by astrophysical backgrounds. 
In preparation for further data from Fermi, as well as new data sets from H.E.S.S., CTA, 
and other future gamma ray observatories, it is worthwhile to explore theories which would 
be expected to result in a prominent line signal. 

Typically, there is some tension between the need to have a large loop level line signal 
without ending up with a large annihilation cross section into continuum photons, which 
could either swamp the line signal, or would over-saturate the thermal cross section. An early 
step in this direction was the "Higgs in Space" model |5j, which was inspired by models of 
dark matter in the setting of warped extra dimensional GUTs (6j[7], and had a Dirac fermion 
dark matter which annihilated through a Z' which coupled preferentially to top quarks. By 
annihilating into top quarks, one could arrange for a thermal relic with masses slightly below 
the top mass which allowed for efficient annihilation in the early universe, but suppressed 
continuum annihilation in cold structures such as galaxies today. 

Ingredients such as s-channel resonant annihilation and forbidden channels in the final 
state have received much attention recently [8], largely in response to the reported excess 
of gamma rays around ~ 130 GeV originating close to the galactic center in the Fermi 
LAT data While interest was inflamed by a subdominant signal at around 115 GeV 

consistent with an additional annihilation channel into jZ 11 13 14 , peculiar features such 
as an excess of photons at the same energy from the Earth's limb 11, 15 16 could point to 
an instrumental origin for the feature. While many of the constructions we consider could 
potentially explain the Fermi signal should it turn out to actually be dark matter, we remain 
agnostic about its origin and consider theories with enhanced lines in general. 

In a previous paper [18] , we considered the forbidden channel mechanism in the context 
of dark matter which annihilates through exotic states with only weak connection to the 
Standard Model itself. We made the important observation that for this class of models, 
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the one-loop continuum annihilation channels cannot typically be ignored, and often play 
an important role in determining both the observability of a line signal as well as the relic 
density. In the current article, we revisit models closer in spirit to the Higgs in Space, where 
the role of the SM-charged mediator is once again played by the top quark. We consider 
UV-complete models, and derive predictions for gamma ray lines, as well as conditions for a 
thermal relic, constraints from direct detection, precision measurements, and the LHC. 
Our article is outlined as follows: in Section [2] we set up the basic module containing the 



necessary ingredients, building on the experience of Ref. 18 . In Section |3j we compute the 
constraints from precision measurements and the observed properties of the newly discovered 
Higgs boson, which help narrow down the viable parameter space, as well as direct constraints 
from LHC searches. In Section |4j we compute the rates for dark matter annihilation at one 
loop, including contributions to the continuum gammas, gamma ray lines, and the relic 
density. Finally, we conclude in Section [6j The appendices contain detailed expressions for 



effective vertices at one loop, complementing those already presented in Ref. 18 



2 Top as a Messenger to Dark Matter 

The basic set-up builds on the models explored in [l8) : the dark matter is a DiracQfermion 
v with no SM gauge interactions, but charged under a U(l)' symmetrjj^J The SM (with the 
subtle exception of the top quark as outlined below) is uncharged under the U(l)', resulting 
in very weak bounds from precision measurements (19). The U(l)' is higgsed by a scalar 
particle, resulting in a massive Z' vector boson and a Higgs-like state $ (which is not to be 
confused the SM Higgs boson h). For simplicity, we assume that any kinetic mixing between 
the U(l)' and the SM hypercharge interaction is small enough to be ignored. Both the Z' 
and $ can act as mediators between the dark matter and some additional fermions ip which 



are charged under the Standard Model. In Ref. 18 , we restricted our attention to the case 
where ip was electroweakly- but not color-charged; this avoided constraints from the LHC 
and prevented annihilation into gluons. In this article, we consider the case where the role 
of ijj is played by an extended top sector. 

We denote by the field which will largely become the right-handed top quark, un- 
charged under the U{1)' and in the (3,1,2/3) representation under the SM (S77(3), SU(2), 
£7(1)) symmetries. Similarly, Q 3 is the usual left-handed quark doublet for the third gen- 
eration. The additional fermions ipL,R have equal charges under U(l)' and are charged as 
(3, 1,2/3) under the SM. By introducing a vector-like pair of ip fields, we avoid disrupting 
the cancellation of any of the purely SM anomalies. The ip fields are also vector-like under 
U(l) ; , and thus no mixed SM-£/(l)' anomalies are generated. The last remaining potential 
gauge anomalies (depending on the charges of v L and Ur) are [/(l)'-gravity 2 or U(l)' 3 , and 
can be arranged by including additional SM singlet fermions which are charged under U{1)' . 
When present, these fermions are essentially irrelevant for the phenomena of interest to us, 
and we will ignore them in the subsequent discussion. 



-'^One could engineer Majorana dark matter, by involving another fermion charged under U{1)' , but this 
would require v-Z' axial interactions, and would generally lead to a large annihilation into gg [18] . 

2 We consider variations where the dark matter v either has vector-like U(l)' interactions, in which case 
one can write down a gauge-invariant mass, or have chiral interactions, in which case the mass for the dark 
matter will have to be generated by the ?7(1 /-breaking VEV. 
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The terms responsible for top-?/' mixing are structurally reminiscent of top-seesaw models 
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Anass = yHQ 3 t R + flll>Ll(>R + Y$ll> L tR (1) 

where H is the SM Higgs doublet, y and Y are dimensionless couplings, and /x is a gauge- 
invariant mass term for if). We define mass eigenstates t and T: 



t R/L \ = ( -sm6 R/L cos9 R / L \ ( i R/L 
Tr/l J V cos °R/l sm6 R/L J \ m R/L 



(2) 



with eigenmasses, 

Mt/r = \[^ + V 2 ^ 2 ) + y2 ( $2 ) T fWS + (/i 2 + 2/ 2 (# 2 > + V 2 ($ 2 » 2 , (3) 
and mixing angles and 0£ given by, 



/x 2 - y 2 (H 2 ) - F 2 ($ 2 > + v/-4/iV(i7 2 ) + (/i 2 + y 2 (# 2 > + F 2 ($ 2 )) 2 
taU0R = WW) 

^ - y 2 (ff 2 ) + r 2 ($ 2 ) + yZ-^v^) + (/x 2 + y 2 (# 2 ) + ^ ($2)) 2 

We identify the lighter mass eigenstate with the top quark of mass m t — 174 GeV, leaving two 
free parameters which we choose to be y(H) and Y(&). Imposing the top mass constraint, 
\x is given by, 

* ~ 1 m\-y\H 2 ) (5j 

and M T is determined We show the contours of //, M T , cos 6r and cos 9l in Fig. [TJ It is 
worth noting that in the limit Y(<$>) — > 0, the mixing angles do not go to zero (unless one 
also requires y{H) — > m t ). For a fixed y(H) ^ m t , the limit Y(<&) — > necessarily implies 
\i — > m t in order to keep the lightest mass eigenstate at m t = 174 GeV. In addition, in the 
SM-like regime where cr, cl <C 1 and the second mass eigenstate is very heavy is obtained 
for y($),/i oo. Finally, the case where y(H) ~ m t and Mt 3> m t corresponds to a value 
of Y(<&) which is parametrically smaller than Mt itself. 

In the mass basis one has diagonal interactions of t and T as well as mixed t-T interactions, 

C = iilp t t + iflpTT+^(c L f-s L i)YW^P L b + h.c (6) 

v2 

+ s L c L \g? R -g? L ]T'fZ lt P L t + h.c 
+ f -fZ'^ gf(s L c L P L + s R c R P R ) t + h.c 

+ ys L SR hit + yc L c R hTT-yh t(c R s L P R + c L s R P L )T + h.c 
- Yc L s R (pit+ Ys L c R (pT T + Y (p t(c L c R P R - s L s R P L )T + h.c 



where 



D tll = dp - i 



(g? R C 2 L + gf L S 2 L )P L + gf R P R ] Z, - igf (c 2 R P R + c 2 P L )^ - ig] R A, (7) 
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Figure 1: Contours of /x, M T , cos6 R and cos9 L in the top partner model. 



and, 



D 



d„ 



(9iA + 9i L cl)P L + gf R P R ] Z, - igl'(s R P R + s 2 L P L )Z^ - igl R A, (8) 



with gj^ = 2e/3, g? = e cos ^ sin ^ (T 3 — (2e/3) sin 2 9 W ) and g = e/ sin9 w the usual 

gauge couplings. In addition, both t and T have the usual diagonal color-triplet couplings to 
gluons, as enforced by SU (3) gauge invariance. We have introduced the short-hand notation 
cl = cos9 L , s R = sin9 R , etc. 



3 Precision and LHC Constraints 

Before turning to the annihilation cross sections and resulting line signals, we examine the 
constraints on the parameter space from precision measurements, LHC searches and direct 
DM experiments. 
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y<H> (GeV) 



Figure 2: (a) Region of the parameter space of Mt and Cl consistent with precision elec- 
troweak data at la (green region) and 95% confidence level (yellow region). Also shown 
is the approximate lower bound on Mt from null LHC searches (dashed line); (b) Contour 
plots for cl,cr and Mt in the ((yH), (Y<&)) plane where the yellow region is excluded by 
EW data. 



3.1 Oblique Corrections 

As a consequence of mixing with T, the Z-t^-th and W-ti-bi EW interactions of the top 
quark are modified with respect to the Standard Model; these deviations are small in the 
unmixed limit of — > 0. Because of the modified EW interactions of top, as well as diagrams 
containing the partner quark T explicitly, there are oblique contributions to the precision 



electroweak observables encapsulated by the Peskin-Takeuchi parameters 24 which provide 



bounds on the mixing parameter cl and Mt- The relevant corrections were computed 



in Refs. 23,25 , where it was found that AS and AU are very small, and the dominant 



constraint is from AT, 

AT = T SM xcl(-(l + S l) + clr + 2sl-^^logr\ (9) 

where, 

r s -r, wi T SM = —^^ c 1.19 . (10) 

For a Higgs mass of = 125.7 ± 0.4 GeV, a global fit to the electroweak data restricts 
AT < 0.14 (0.10) for AS ~ AU ~ at the 95% (68%) confidence level (26) . For a fixed 
value of the mixing cl, this corresponds to an upper bound on Mt- The constrained region 
is plotted in Fig. |2Fa), and indicates that for cl ~ 0.1, masses Mt < 3 TeV are consistent 
with precision measurements. Also shown on the figure is the line corresponding to Mt > 
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Figure 3: ft, and 99 production by gluon fusion and their decay into 77. 
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Figure 4: Contours for Y^gg/Y^ . The yellow region is excluded by precision electroweak 
measurements. 



475 GeV, which is a conservative estimate for the lower bound on the mass of the T quark 
based on a collider searches (see below). 

We show in Fig. |2^b) the contour plots for cl, cr and My in the parameter space of y(H) 
and Y with the region excluded by precision measurements overlaid. In the allowed 
region, cl is smaller than ~ 0.2, but cr may still be sizable provided T is not too heavy. 



3.2 Constraints from Higgs Measurements 

In addition, there can be constraints on cl and cr from production of the newly discovered 



Higgs boson 29 at the LHC, which receives modifications due to the shift in the top cou- 



pling to the Higgs, as well as corrections from the T quark itself, which contributes to Higgs 
interactions with gluons and photons (see Fig. |3j. However, in the region consistent with 
precision electroweak constraints, Higgs phenomenology remains almost very close to Stan- 
dard Model-like at the LHC. For example, the corrections to h — > gg (which also characterize 
the shift in inclusive Higgs production) are, 



SLSRFi/ 2 (r t ) c l c r F 1/2 (t t ) 
m t /v ttit/v 



oc 



Fi /2 (r t ) 



m t /v 



where the function Fi/ 2 of 



m\/ (4m^) can be found in Ref. 



30 



As illustrated in 



Fig. [1J these are much less than 1% deviations for the parameter space of interest. Such 
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tiny deviations might be measurable at a future linear collider and/or Higgs factory, but are 
unlikely to be accessible at the LHC. Similarly, the correction to h — > 77 is likely to be even 
smaller since it is dominantly mediated via a W loop, and would be extremely difficult to 
detect. 

One could also search for the $ boson, which, once produced, would decay via its mixing 
with the Higgs. However, the t and T loops contributing to $ production at the LHC almost 
cancel since tan 9 a/ tan#£, = m t /m T , and 



^ip^gg OC 



CLSRFi/ 2 {r t ) s l c r F 1/2 (t t ) 



m t /v 



2 2 

C L CRS R r 

oc \Fi/ 2 {T t )\ 

m t 



1 



Fi/ 2 (tt) 



F 1/2 (n) 



resulting in the production rate of $ being very suppressed. 



(12) 



3.3 Constraints from LHC Searches 

The primary constraints from the LHC are from the null searches for exotic quarks such as 
our T. As a vector-like quark which mixes with the top, the primary decays are expected 
to be: T — > Wb, T — > Zt, and T — > ht. In addition, if either $ or Z' are light enough, 
one could have decays into either one of them plus a top quark. Null searches for the Wb, 
Zt, and ht channels at the LHC result in a conservative (in the sense that larger masses are 
permitted for any configuration of branching ratios) bound of My > 475 GeV 27, 28 . In 



the case where T decays entirely into th, the current bound derived from ATLAS data is 



425 GeV 27 



In addition, there are constraints on various types of Z' bosons that may be inferred 
from LHC searches. The primary interaction of the Z' with the Standard Model is with 
top quarks. As we shall see, the result is suppressed production cross sections which are 
generically safe from LHC bounds. 

Given the large coupling to top quarks, the dominant Z' production is via gluon fusion 
through a top loop, as in Figure |5j Detailed expressions for the effective Z'-g-g vertex are 
provided in Appendix [Bj The Z' can subsequently decay into vv leading to a missing energy 
signature, or back into top quarks leading to a resonance in the invariant mass distribution 



of tt pairs. An interesting feature of this signal arises because the Landau- Yang theorem 31 
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insures that the amplitude must vanish when the Z' is on-shell. As a result, the signal is 
much smaller than one might naively expect, and would have an interesting shape around 
the Z' resonance if it were visible and the Z 1 width were large enough to permit it to be 
experimentally mapped out. However, from the point of view of constraining the parameters, 
one can derive more stringent bounds by considering the process where one of the initial 
gluons is highly off-shell, resulting in a final state consisting of the Z' and an additional 
hadronic jet (and after the Z' decays to ti leads to a tij signature). 

The choice of px cut on the associated jet can be optimized based on the need to pull 
one of the initial gluons off-shell so that the Z' production rate escapes from the Landau- 
Yang suppression, while at the same time not choosing a cut so stringent as to result in a 
large suppression from the parton distribution functions. Fully optimizing the search cuts is 
beyond the scope of this work, but for a choice of pr > 30 GeV, large enough to suppress the 
ti background and such that collinear logarithms do not spoil perturbation theory, we find 
that even for large couplings, the deviation at yfs = 8 TeV is never more than about 0.1% 
of the SM tij rate^J with the largest potential deviations occurring when the the Z' mass 
is just above the ti decay threshold. For larger values of the px cut, the deviation is even 
smaller. Such deviations are well within the uncertainties in the SM rate from theoretical 
uncertainties such as the gluon parton flux and scale variation, and thus will be very difficult 
to extract from the background. This situation is essentially unchanged at = 14 TeV. 

When the Z' decays directly into dark matter, it produces a dark matter mono-jet sig- 
nature 33 . For very heavy Z's, well-represented by a contact interaction, the current LHC 
bounds are a factor of few too weak to bound our construction [34], but future prospects 
in this channel would appear to be good. As the Z' gets lighter and is no longer well- 



approximated by a contact interaction, these bounds are typically somewhat relaxed 35 



The Z' could also be produced through its induced electroweak interactions, through 

in 



processes such as qq — > V* — > Z'V (where V, V = 7, Z), which was considered in Ref. 36 
the context of the Fermi 130 GeV line feature. For Z'-Z-j couplings large enough to explain 
the Fermi feature, the LHC can provide useful bounds (though it is somewhat questionable 
whether such large effective interactions are consistent with a perturbative UV complete 
theory) in the regime of weak DM-Z' couplings, but in our construction the effective vertex 
is small enough that there is essentially no useful bound from this channel, even for very 
light Z's. 

A Z' with large coupling to top quarks can also be produced as radiation from a ti pair 
through processes such as gg —> tiZ' . Its subsequent decay into dark matter leads to a top 



pair plus missing energy signature 37 , or it can decay back into top quarks leading to a 
four top final state 38 -41 . For either process, the rate is modest because of the phase space 
suppression and is unlikely to be a serious constraint from current LHC data. However, 
prospects are good for observation at higher energy LHC running. 



3.4 Constraints from Direct Detection 

For low Z' masses, constraints from null results of direct searches for dark matter (e.g. from 
Xenon-100 |42|) can be potentially rather strong. In particular, it implies a bound on the 

3 When the Z' can decay both into vv and ti, the branching fraction for Z' —> ti is less than 10 -2 , since 
the top coupling is suppressed by c 2 R compared to the v coupling. 
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Figure 6: Cross section for annihilation into the 3-body final state tWb for several values of the 
relative velocity v and Mz> mass, as indicated. We have set cl = cr = 0.2 for simplicity. 



degree of kinetic mixing between the Z and Z' , since that parameter controls the coupling 
to light quarks^} As the direct scattering of v with nuclei is entirely on protons rather than 
neutrons, the Xenon-100 bounds (which assume equal scattering with both protons and 
neutrons) require a pu < 3 x 10~ 44 cm 2 , corresponding to the restriction 77 x g z > < 1CT 3 [5], 
where r] is the coefficient of the FyF' term in the effective action. 

In principle, rj is a free parameter which is additively renormalized by loops of the t and 
T quarks. If the U(l)' is embedded in a non-Abelian symmetry, it will naturally vanish at 
the scale at which that symmetry breaks down to U(l)' (which acts as a UV-cut-off on the 
log divergence present in the current framework). Its natural value thus depends on both 
the Z 1 gauge coupling and the symmetry-breaking scale. For the parameters of interest to 
us, we find r)/gz< < O(10~ 3 ), and thus consistency with the direct bounds requires at most 
very modest tuning of rj. 

3.5 Flavor Considerations 

Finally, we should mention the possibility of flavor violation. We have chosen to identify the 
SM quark with which the messengers ipL,R mix as the top quark, but there are no symmetries 
that prevent similar mixing terms with the up and charm quarks. Indeed, the CKM rotations 
would almost certainly induce such terms after electroweak symmetry-breaking. From a 
purely phenomenological stand-point such mixing must be small in order to avoid constraints 
from D meson mixing and rare top decays such as t — > cZ'* — » ccc. We will invoke parameters 
such that these decays are unobservably small, but they represent a generic feature and an 
interesting signature that is worthy of further investigation. 

There is also direct scattering through the induced coupling to gluons mentioned above, but it leads to 
predictions well below the current bounds. 
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4 Dark Matter Annihilation 



4.1 Relic Density 

In the forbidden channel scenario, the correct relic density can be obtained for DM masses 
slightly below the mass of the particle running in the loop. For our case this particle is the top 
quark, favoring dark matter masses between 150 - 170 GeV. While such masses will typically 
have the correct ballpark relic density, particularly when the annihilation occurs somewhat 
close to the Z' resonance (~ 350 GeV), the detailed story depends quite sensitively on the 
entire set of model parameters. For example, annihilation into the 3 body final state tWb 
can be important (see Fig. [6]), especially when the dark matter has sufficient velocity to allow 
the Wb system to form a close to on-shell top quark. For regions of parameters, processes 
at the loop level can also be significant. For example, as we will see below, annihilation into 
Zh is often the most important channel both today and in the early Universe. 

In addition, there are viable points for which the Z' is only slightly more heavy than the 
dark matter and the relic abundance results from annihilation into Z'Z'. For example, one 
has the correct relic abundance for M' z ~ 120 GeV when M v = 100 GeV. Note that the 



bound on 77 imposed by direct detection (described in section 3.4) makes the contribution 
from annihilation via Z-Z' mixing irrelevant for the relic abundance calculation in virtually 
all of the interesting regions of parameter space. 

4.2 One Loop Annihilations and Gamma Ray Lines 

At one loop, we have a plethora of final states available as shown in Figures [7] and |SJ In 
addition to line signals from jZ, ^Z' and 7/1, there is potentially an extra line due to the 
7$ channel. Moreover, the one-loop annihilations Zh, Z'h, Z§, WW, ZZ, ZZ' and bb 
can contribute to a sizable continuum of gamma rays. The analytic results for the effective 



vertices can be found in the appendices of Ref. 18 and (for annihilation into bb) appendix 
|A} Annihilation through the Z'-^-h, and Z'-h-h effective vertices vanishes. 

Compared to the original Higgs in Space model [H], there are significantly more one 
loop diagrams. In that model, 1 and only t R couples to Z' while the T-b-W, t-T-Z, 

h-T-T and &-t-t couplings are suppressed and m T is large. In the current, UV-complete 
construction, cl and cr are related and go to zero together. In the limit of small mixing 
between tR and ip, the t-t-Z' coupling is suppressed and while T has large couplings to both 



Z and Z' they are both largely vector-like and thus do not contribute 18 . However, there is 
a new feature in the present work that arises because of the "mixed loops" where both t and 
T run in the same diagram. The mixed Z'-t-T couplings are less suppressed for small mixing 
angles and in addition there are enhancements due to the presence of two different mass 
scales in the loop, leading to a sizable annihilation rate into the vector + a scalar channels, 
even for modest mixing angles. 



To avoid the overwhelming 1-loop annihilations into gluons 37 43 44 , we focus on the 
case where DM has vector-like couplings with Z' . Although we are forced to live in a regime 
of relatively small mixing angles by the EW precision constraints, we can still have large 



annihilation cross sections into gamma ray lines, as shown in Figures [9 10 and 11 



In Ref. 18 , where no fermion mixing was involved, we found that the dominant channels 
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Figure 7: One-loop Z' effective vertices induced by t and T couplings to 7, Z, Z', /i, 
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were Z<3> and 7<Hj For the top quark mediator model, we now have the possibility to produce 
the SM Higgs in the final state. Moreover, because of the presence of "mixed" loops involving 
both t and T, we find that Zh dominates over 7/^] and the discrepancy grows in the limit 
of small mixing where g a z' *C gtrz 1 - 111 most of the (y(H), Y (<$>)) parameter space, Zh is 
the dominant annihilation channel, although as shown in Figure 10, there are regions of 



parameter space where the Zh continuum is relatively suppressed, resulting in a large 7Z 
over continuum ratio. 

Roughly speaking, we can classify the situation as to whether the dark matter is heavier or 
lighter than [Mz + Mh)/2. For the heavier case, the dominant one-loop annihilation channel 
is Zh. As discussed above, the relic abundance will typically work out for dark matter 
masses somewhat below the top mass, where there is significant phase space loading of the 
Zh annihilation rate. Still, a reasonably large line to continuum ratio, as large as O(0. 1) 
can result, leading to striking line features as illustrated in Fig. 11 In the lighter case, the 



one-loop continuum tends to be dominated by WW (which is always significantly larger 
than ZZ due to two very different masses mt vs. in the loop, leading to some amplitude 



J This result can be understood as follows: In minimal anomaly-free constructions realizing axial Z' 
couplings, 1-loop annihilations into two gauge bosons are suppressed due to a cancellation between 1-loop 
diagrams. Such cancellation does not occur for the 7$ and Z$ channels, which end up dominating. 

6 The "mixed" loops are larger than the top loop, but for the 7/1 channel, they cancel each other exactly 
as t and T have equal couplings to the photon. Therefore, the 7/1 cross section is entirely controlled by the 
t loop. On the other hand, for Zh annihilation, the 'mixed" loops do not cancel each other, as gtz 7^ grz, 
and dominate. 
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Figure 9: Benchmark predictions for cross sections (as labelled in the plot) and the ratio of line 
to continuum rates, as a function of the Z' mass. Note that the predictions corresponding to the 
center panel are essentially unchanged if M v = 150 GeV rather than M v = 135 GeV, for the same 
choice of (yH) and values. 
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Figure 10: Predictions for line signal rates and ratio to continuum for two choices of DM masses 
in the plane of (y(H),Y ($)). The shaded region is excluded by EW precision tests. In the bottom 
plots, the choice of the DM and Z' masses guarantees that the correct relic abundance is obtained 
for essentially the whole plane. 

enhancement) which is nonetheless is subdominant to the line annihilation into ^Z, resulting 
in a line to continuum ratio of order one (see Fig. [9]). The correct relic abundance can be 
obtained from annihilation into Z'Z' with M' z ~ 120 GeV if M v = 100 GeV. 

Finally, note that a line from annihilation into 7Z 7 is possible, although only for a rela- 
tively light Z' which will not typically lead to the correct relic abundance. In addition, there 
will typically be a very large one- loop continuum contributions from Z'h and/or Z$. The 
7$ line may be as intense as the jZ one, however, there is typically a relatively small Z' 
mass window where 7$ is kinematically open while Z<3> is not. When Z$ is open, it tends 
to dominate all channels by several orders of magnitude. 

To summarize this section, we find that when the top quark acts as the primary mediator 
between dark matter and the Standard Model via an s-channel resonance, the dominant line 
signal tends to be 7Z with 7/1 relatively suppressed. The one loop contributions to the 
continuum annihilation cannot be ignored, and can be dangerously large from annihilation 
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Figure 11: Example of a gamma-ray spectrum obtained for (yH) = 177 GeV, {Y3>) = 200 GeV, 
9z' = oiz 1 = ^' leading to a two-line signal at = 128 GeV from 7Z annihilation with ajz v = 
2 x 10~ 27 em 3 s _1 and at E 1 = 71 GeV from 7$ annihilation with a^v = 1.35 x 10~ 28 cm 3 s -1 . 
The dominant continuum is from WW and Zh annihilations, ay/w v = 5.2 x 10~ 27 cm 3 s _1 and 
azh,v = 1.7 x 10~ 27 cm 3 s -1 . 



into gg if the dark matter has axial vector interactions, and completely swamp any line signal 
if Z$ is open. 



5 An Alternate UV Completion 

Large line signals require axial couplings in the loop [18]. The current construction realizes 
this need by mixing with the SM top quark, making use of its chiral nature, but its effec- 
tiveness is limited by the constraints from precision data, which require at most a modest 
level of mixing. In this Section, we consider a simple extension of the basic UV completion 
presented above in which the heavy partner fermions ip have axial couplings with the Z' 
even before mixing. Consequently, there are large line signals even in the "no- mixing" limit, 
where the top coupling to the Z' coupling goes to zero (preventing annihilation into tt), and 
allowing for dark matter masses well above the top mass. 

We start with the previous UV completion of Section [2] and include a second vector-like 
doublet ip 2 = (^2L,^2i?) = {^l^r) in addition to the original ipi = (Vhi,, 'VW) = {^l^r)- 
Both carry the same SM gauge charges as t R , and ipi also carries U(l)' charge q\ — g$ = 1, 
whereas ip2 has q2 = 0. For simplicity, we set the vector-like masses of ip\ and 1/^2 to zero 
and assume that the Higgs Yukawa coupling between Q 3L and ip 2 R is negligible. We have 
Yukawa interactions: 

C = y HQ 3 t R + Y ^lt R + A ^F L r R + A' (13) 

We recycle the formulae Eq. for the mass eigenstates t and T, with the replacement 
that the \x parameter is replaced by A($) in the formula for tan6^/£. Imposing that the 
lightest mass eigenstate has the SM top mass fixes the value of \i = A($). In addition to 
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t and T, there is another massive state T" = ip + that does not mix with tn and which has 
mass given by Mt> = A'($). In this way, we realize the existence of two massive states T and 
T' with axial-vector couplings to the Z' , even in the limit that the SM top quark is purely 
unmixed. 

With respect to the previous case, the couplings to the Z 1 and $ are modified, whereas 
the couplings to the Z and SM Higgs are not. In particular, t# no longer couples to the Z'. 
The modified couplings can be summarized: 

V 9Z> 2 A 9Z> 2 V 9Z' 2 A 9Z> 2 / 1A \ 

9ttz' — 2 ° L) ^ ttz ' ~ 2 Cl ' ■S'ttz' — 2 Sl ' S'ttz' _ 2 Si ' 

£trz' = "2~ CiSi ' &rz' = — 2~ ClSL ' ^ ^ 

2/^ = -Yc L s R + \c L c R , y^T^ = Ys L c R + Xs L s R , yf m = 0, y%r$ = (16) 

(17) 



= ^[^(clC/? - s L s fi ) + A(c L s i? + s L c R )], y t y # = ^[^(clc^ + s L s R ) + A(c L s fl - slc^)] 



and in addition 

n v 9z ' n A - 9z ' 

9t't'z' — ~^-> 9t't'z' — ~^-> \ 18 ) 

Vt'T'® = A , Ut'T'® = (19) 

To simplify the parameter space, we set all pseudo scalar couplings of ip to $ in the 
Lagrangian to zero, although the mass eigenstates end up with pseudo-scalar couplings 
for the $-t-T vertex. Having set the DM pseudo-scalar couplings to zero, the one-loop 
diagrams mediated by $ will be p-wave suppressed, leaving the relevant contributions from 
Z' exchange. 

As mentioned above, this framework evades large annihilation into tt when the dark 
matter mass is above the top mass. However, because the constraints from precision mea- 
surements are greatly relaxed, T and T" can as light as the direct search bounds of ~ 475 GeV, 
and one can obtain the correct relic abundance from annihilation into tT/Ti, TT, T'T 1 or 
One-loop DM annihilations are dominated by the contributions of T and T' , although 
mixed (t,T) loops still play an important role. Zh dominates the continuum when it is 
kinematically accessible. 

There are line signals from jZ, jZ' and 7$, whose relative sizes depend on the choice 
of parameters, in particular the relative masses of v, Z' and $. We illustrate some cross 



section predictions in Figs. [12] and [L3| Typically, we find that the dominant line is due to 
7$. In the scan plots of Fig. 13, we therefore show contours for <T 7 $f and the ratio cr 7 $t> over 
continuum, finding that the value of cr 7 $i> is essentially constant in the plane (yH), (Y<&), 
as it is nearly independent of the mixing angles. However, the other channels depend very 
sensitively on them. As a result, the line over continuum ratios range over several orders of 



magnitude depending on where in the plane one considers. In Fig. 14 , we show an illustrative 



2-peak spectrum due to M v = 500 GeV DM annihilation into 7$ and 7Z 7 . 
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Figure 12: Predictions for cross sections and line to continuum ratios in the second UV completion. 
The upper plot corresponds to a light DM case where the 7$ line signal is huge and the continuum 
is essentially absent for Z' masses leading to the correct relic abundance. The center plots are quite 
typical of what happens for a generic choice of parameter, with Zh and Z§ channels dominating. 
The lower plots correspond to the case where (yH) —> mt and mixing angles vanish. In this situation 



ZQ and 7$ are the main channels 18 . 
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Figure 13: Predictions in the second UV completion, scanning over the plane of y{h)-Y{&), where 
the shaded region is excluded by EW precision tests. In the upper plots, the choice of the DM 
and Z' masses guarantees that the correct relic abundance is recovered in the whole plane through 
vv —7- Z 1 Z 1 . For the lower plots, the relic abundance is controlled by annihilation into T'T' (M?* ~ 
530 GeV), Z'§ or tT in the region Mt < 800 GeV, in the early Universe, with the correct relic 
abundance recovered in a significant fraction of the ((yH), (¥$>)) plane, (particularly in the region 
where the line over continuum ratio is larger than 10~ 2 ). 

6 Summary 

We have considered a class of models where the dark matter is a Dirac fermion whose con- 
nection with the Standard Model is primarily through the Top quark. We considered two 
variations of UV complete models, and found that both lead to large gamma-ray line sig- 
natures from annihilation into either ^Z or 7$, where $ is the scalar responsible for the 
breaking of a new U(l)' gauge symmetry. In the first variant, the Z' has sizable coupling to 
the top quark and a mass only slightly below it, realizing the forbidden channel mechanism 
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Figure 14: Example of a gamma-ray spectrum obtained in the second UV completion for (yH) = 
174.6 GeV, Y = 0.1, A' = 0.62, g z > = 1 = g^ z ,/2, leading to a two-line signal at E y = 164 GeV from 
jZ' annihilation with a^z'V = 4 x 10~ 29 cm 3 s _1 and at E~ f = 463 GeV from 7$ annihilation with 
a 7 $v = 2 x 10~ 28 cm 3 s -1 . The dominant continuum is from Z& annihilations, a z $v = 5 x 10~ 28 
cm 3 g -i Yi e n c abundance is controlled by vv — > T'T' which is open in the early universe but not 
today. 



by which the dark matter annihilates primarily into top quarks in the early universe, but this 
annihilation is shut off for cold dark matter particles today. In the second UV completion, 
the Z' coupling to top quarks is very suppressed, and the dark matter mass can be much 
larger than the top mass. In both cases, Zh (and, depending on the $ mass, Z$) is typi- 
cally the dominant annihilation channel, illustrating the need to carefully consider one-loop 
contributions to the continuum photons in theories with enhanced line signals. 

Gamma ray line signatures are an important item on the menu of indirect searches for 
dark matter. Models such as the ones considered here, which suppress the continuum from 
annihilations today, are particularly well-probed by these searches. What is perhaps more 
surprising is the fact that these models are surprisingly immune to constraints from the 
LHC and direct detection as well, which would suggest that the gamma ray sky hides unique 
probes of interesting dark theories. 
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Figure 15: Topology for the effective vertices. 



A One-loop annihilation into bb 

In this appendix, we summarize the one-loop expressions for the effective vertex and amplitudes- 
squared needed for the bb channel. One-loop expressions for all other channels are listed in 
the appendix of Ref. 18 . The topology for the loop diagrams considered here is shown in 
Fig. 15 We express the amplitude in terms of two-point (B ) and three-point (Co) scalar 



integrals where: 



Co = C (MlMl4M*-ml,ml,m 2 3 ) 

£? (23) = B^Ml-mlml), 

S (13) = B {4M*;ml,ml), 

£? (12) = Bo{Ml-m\,m\) . 



(20) 
(21) 
(22) 
(23) 



In the following expressions, mi = 772,3 = m f an d rri2 = My/- However, to derive our results, 
we have computed all of the expressions in these appendices in terms of the general masses 



as depicted in Fig. 15 



The effective vertex for Z' a — > bb is given by the expression: 



V 7.1- 



(24) 
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where the loop coefficient A is: 
9 



2 



A = 32^4 ^N^Ki - 4) - ™? (^4 + ^»(^o(12) + B (13) 
-2S (23)) + 8C M*{vf f - of}) J - 2M^ Uaf} - t$) (s (12) + B„(13) 
-2^£?o(23) + 8C M^ + C mJ(5af} - 3vf f )^j + M^^m 2 f (3af f - - vfj) \B Q {12) 
+S (13) - 2 ^ (23) + 8C M^ ^ + 8M 2 (afj - vfj) \-2B {l2) - S (13) + 35 (23) 
+8C M 2 + 1 J + 8C Q affmj J + 4M>J(af}(25 (12) + 5 (23) - 1) + vf f (-2B Q (12) 



+B (23) - 1)) - 2C m f (afj + vfj) + 4C M^(af} - vfj^j . 



(25) 



The matrix-element-squared takes the form: 



lMl = Tj^j* { 4MP4 + W I (26) 

B t£ production mediated by Z 7 

In this appendix, we provide the expressions for Z' production at the LHC The Z' is 
produced via loop coupling to gluon pairs through the top sector. Similar to the ^'77 



process, the Z'gg vertex convention is shown in Fig. 16 and is given by [18 



Kz> = E i 1 + 2 ™^o(Mi„ 0, 0, m}, m% m})) , (27) 

/ 

where N[ is the color factor of the internal fermion, g s is the SU(3) coupling, ajj is the 
axial-vector Z' — f — f coupling. The tensor structure of the ggZ' vertex is given by 



2~M\ 

The loop function Cq can be reduced in this instance to 



C (M 2 z „0,0,m 2 ,m},m 2 f ) = ^ log 2 {^-^ ) , (29) 

4m 2 , 



where /3f = \[1 — j^f. Thus, the effective vertex becomes 



v ^a _ r ^a § ab^ N £9saff ( \ m ) , 2 ( ' Pf ~ 1 \\ ( o Q) 

V ggZ .-l 9gZ ,d 2^ 4?r 2 { l + Ml i0g \p f + l))' m 
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Figure 16: Effective Z'gg vertex mediated through t and T loops. 



The Landau- Yang theorem forces a vanishing vertex when any of the external bosons 
is on-shell. Therefore, the leading production mechanism of the Z' is in association with a 
gluon or quark shown in Fig. [5| The Z' subsequently decays to it pairs. 

The amplitude for g(qi) + gfa) —> t(pt) + i(p' t ) can be expressed as: 



(-9aP + 



Pz', a Pz',p 



AW* = ^ (4 + 4' 75 ) v ( P ' t) {m ; ri ;_ M| X r / r;g tWttfe ), (si) 

where we have expressed the Z'gg vertex as Vz> gg = AT^, and e(?i)^ and e(q2) v are the 
polarization tensors of the incoming gluons. Note that when we square this amplitude and 
average over the gluon polarizations e(gi)^ and e(q2) v , we must do so in such a way that 
only the physical (transverse) polarization states contribute to the matrix-element-squared. 
To do this we adopt the general prescription: 

£ e„( ft) AOvfe, Ai) = + + - , (32) 

where i = 1,2 and the arbitrary vectors rij have to satisfy the relations: 

^(ft, Ai)e*,(^,Ai) = , nf £ ^(ft, Aj)e*,(ft, Aj) = (33) 



A 1 

n7 



A=l,2 A=l,2 
,2 



together with 7^ and n\ 7^ n 2 . We choose ni = q<i and n 2 = qi such that: 

A<=1,2 



£ ^vCft, A,) = + 2 ^v + ^v ; (34) 



where s = {q\ + g 2 ) 2 - 

Finally, summing over the final-state top quark polarizations and averaging over the 
initial state gluon polarizations, we find for the matrix-element-squared: 



Y)M\ 2 - ( a a) 2m t\ A \ 2s 



32 (( s - miy + Mf,r|,; 



Ml 



Z'J 



(35) 



We see that in the limit where the Z' goes on-shell (s — > M|,) the amplitude-squared vanishes 
in accordance to the Landau- Yang theorem. 
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